Production of bright entangled photons from moving optical boundaries 
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We discuss a new mechanism of generating two separable beams of light having a high degree of 
entanglement in momentum using a fast and sharp optical boundary. We identify three regimes of 
light generation depending on the number of resonant interactions between the optical perturbation 
and the electromagnetic field. The intensity of the process is discussed in terms of the relevant 
physical parameters: variation of refractive index and apparent velocity of the optical boundary. 
Our results also open the possibility of generating entangled light by different means of exciting 
zero point fluctuations such as the dynamical Casimir effect, and since it results from a two-mode 
squeezing it is very robust against temperature on a wide range of wave lengths. 

PACS numbers: 



Many of the theoretical schemes and experimental ap- 
plications being proposed and developed in the context of 
Quantum Information (including quantum computation 
and information processing teleportation [2j, etc.) 
rely on the possibility of generating entanglement be- 
tween different quantum systems. Though entanglement 
can arise in nature even from the simplest interactions 
and even at high temperature [!, the degree of en- 
tanglement achieved is usually very small. Photons are 
probably one of the most robust quantum systems where 
entanglement can be created and manipulated || , hence 
they are among the best candidates for future practical 
applications. Presently, entangled photons are produced 
experimentally via parametric down conversion (PDC), 
which is a nonlinear process having small efficiency [6|]. 
Hence, it is useful to identify new mechanisms of gen- 
erating entangled light and to develop sources of entan- 
gled photons with high brightness to improve contrast in 
present day and future experiments 0]. 

Recently, a series of papers have introduced the 

concept of Time Refraction (TR), which describes how 
the classical and quantum properties of light are al- 
tered by the sudden change of the optical properties of a 
medium. TR results from the symmetry between space 
and time, extending the usual concept of refraction into 
the time domain. Like the Unruh effect [icj |. the Hawk- 
ing mechanism 11 1 and the dynamical Casimir effect 
the quantum theory of TR predicts the excitation of vir- 
tual particle from the turmoil of Zero-Point Fluctuations 
(ZPF) and the emission of pairs of real photons, each 
emitted in a different direction. The number of pairs 
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emitted is proportional to the variation of the refrac- 
tive index of light associated with the optical pertur- 
bation. For any realistic experimental parameters, the 
mean photon number produced in the optical domain 
from the vacuum state is smaller than 1 and therefore 
a better approach was needed to achieve a measurable 
effect. In a recent work [13|, a different process of excita- 
tion of ZPF was proposed using a non-accelerated opti- 
cal boundary moving with apparent superluminal velocity 
across an optical medium. Like TR, this effect also leads 
to the emission of photons pairs, but now the moving 
optical boundary works as a relativistic partial mirror, 
producing a considerable Doppler shift, altering radically 
the intensity of the interaction between light and matter, 
yielding a potentially measurable number of photons by 
choosing adequately the velocity of the optical bound- 
ary. The pairs of photons predicted in [f| and [HI are 
highly correlated: a photon cannot be emitted without 
the emission of its counter part. This paper addresses 
the emission of radiation produced by a non-accelerated 
sharp optical perturbation, moving with apparent super- 
luminal velocity, extending the results in [13[ from a one 
dimensional to a three dimensional geometry and analyz- 
ing the entanglement of the photon pairs emitted. 

We begin by considering an infinite optical medium 
which is swept by an optical perturbation, described as 
a sharp variation of the refractive index of the medium 
with apparent velocity u (see Figure 1). Optical per- 
turbations with superluminous profiles can be produced 
via the interaction of a short laser pulse with a plasma 
14 1 or by a laser or electron beam sweeping across a gas 
15j . In this context, the apparent velocity u describes a 
delay of the change of refractive index between different 
points of space and does not refer to an actual velocity of 
propagation of the optical profile. Henceforth u can take 



2 



values arbitrarily large, even larger than c. 

Following the method developed in jl3j |. we describe 
the process of interaction between the ZPF and the op- 
tical perturbation in a reference frame S' , with velocity 
Voc = —c 2 /u < c relative to the laboratory reference 
frame S, where this optical boundary will be perceived 
as moving with a infinite velocity: 
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As a consequence of the relativistic phase invariance, the 
refractive index of the medium in the S frame and in the 
S' frame (respectively n and n') are different [§]: 
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where /3 = v^jc — —c/u, 7 = (l — 2 ) and 6 is the 
angle between the velocity of the optical perturbation u 
and the wave vector k . 

In the S' frame the problem is identical to a TR and 
can be solved by imposing the continuity of the dielec- 
tric displacement and the magnetic induction fields and 
corresponding field operators during the time dis- 
continuity of the refractive index, or equivalently, by im- 
posing phase matching conditions at the optical bound- 
ary. Back in the S frame, the optical perturbation can 
be perceived as a four port device, coupling two initial 

complex plane wave modes: (fii(r) = exp[— i k j. r ] and 
(j)a(r) — exp[— i k a .~r*] existing for r > ut, with two 
final complex plane wave modes <fit(r) — exp[— i k t- r ] 
and r (~r > ) = exp[— i k r . r] existing for ~r* < ut, which 
satisfy: 
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-hi 2 [f + 9t] ~u 1 1 - hsm0 l ~u±, (5) 
where / = (cos^ - P/th), A it = hi/h t , A ir = hi/h r , 

1/2 

hi,t, r = [7 2 (cos6>j - P/n it t, r ) 2 + sin 2 9 t ] , g t>T = 0(1 - 
/3nt, r cos di)/rit t ri 6% now is the angle between the veloc- 
ity of the optical perturbation u and the wave vector 

fcj. The different values of the refractive index m, n t , 
n r and n a for the incident, transmitted, reflected and 
anti-incident waves take into account the possible dis- 
persion of the optical medium prior and after the optical 
perturbation has passed. 

Like equation ©, equations (O to © are also derived 
from the invariance of the phase of light between any two 
different inertial frames [9| and correspond to a double 
Doppler shift. For values 9i ^ 0,7r, (fit and <fi r are calcu- 
lated as: 

9 t = arctan [sin^/7 2 (/ + A it g t )] , (6) 
r = arctan [sin 9i/-f 2 (/ — A ir g r )] . (7) 




Figure 1: Space-time schematic diagram of superluminal 
space-time refraction. In the S frame (bold), the optical per- 
turbation is observed as moving along the x axis (from left 
to right) with apparent velocity u whereas, in the S' frame 
(dashed), the optical perturbation alters the refractive index 
of the medium simultaneously for all point of space at instant 
t' = 0. Both frames can be related using a standard Lorentz 
boost with [3 = Voo/c = —c/u < 1. 



These expressions correspond to the generalized Fresnel 
formula for a moving superluminal partial mirror. Notice 
that \9t\ is mainly limited to a solid angle around the 
direction of u (see Figure 2). 

Using the continuity conditions for dielectric displace- 
ment and the magnetic induction fields at time 
r = tu , the annihilation and creation operators for 
these modes can be related as: 
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where A = (l + a 2 ) /2a, B = (l - 1 

1 /2 

a = [niG(0i,ni)hi/n t G(0t,rit)ht] and G(9,n) 
(3ncos(6), satisfying A 2 - B 2 = 1. 

As demonstrated in references QjJilli], the two- mode 
squeezing transformations (JSJ and @ imply that, after 
the optical perturbation has passed, an initial vacuum 
state can be expressed in terms of the new eigenstates of 
the field according to: 



|0),|0) Q = £ n C7„ 
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with C„ = 1 1 - \z\ ) z n and z = B/A. Equation (JTT 
implies the emission of pairs of photons moving along 
different directions with wave vectors k t and k r , ac- 
cording to equations © , ^ and © . The mean number 
of photons with wave number k t and k r is : 
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Figure 2: Emission spectrum: Angular distribution of the 
number of photons emitted by the optical boundary corre- 
sponding to a change of refractive index from nt — n a — 1.1 
to nt — n r = 1.5 and for /3 = 0.9 (dashed) and for f3 = 0.99 
(bold). Notice that the emission predominantly concentrated 
in a very small solid angle around the resonance angles. With 
jfl = 0.9 there is only the resonance for /3riiCOS#[ es — > 1, 
whereas with /3 = 0.99, both resonances /3n, cos B\ ea — > 1 and 
/3n r cos 9r es — ► lexist. 



clearly identify B T t es and #£ es , calculated from #[ es using 
equations © and respectively. 

Equations © and © change the vacuum state into 
another pure state for which the entanglement entropy is 
the canonical entanglement measure [25j. The entangle- 
ment entropy Evn{pss') of a system composed by two 
sub-systems S and S' and described by a density operator 
pss' is calculated by tracing out one of the sub-systems 
in the density operator, i.e. p s = Tr$> [pss'] and calcu- 
lating Evn{psS') = ~Tr [p s In p s ~\ . In the present prob- 
lem the subsystems are the two distinct sets of photons 
emitted corresponding to <f> r and <f>t, yielding 

E VN - In (1 + (N t )) (1 + {N t )) - (N t ) In (N t ) . 

Notice that Ey n is basically the Shannon entropy intro- 
duced by increasing a photon pair in the system. For 
z — ► 1, the entanglement diverges as the system ap- 
proaches the resonance condition and the maximal en- 
tanglement state is achieved, i.e. 
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According to Eq. (fTTj) the number of photons emitted 
diverges for z — > 1. In the one-dimensional case studied 
in [13[ this could only be achieved if a perfect match- 
ing between the velocity of the optical perturbation u 
and Hi such that /3rii = crii/u — > 1. However, in the 
three-dimensional case there is an extra degree of free- 
dom corresponding to the angle between the wave vector 
and the direction of the apparent motion of the optical 
perturbation, and the condition z — > 1 can be achieved 
for both (3rii cos Q\ es — > 1 and/or f3n r cos r r es — > 1. Un- 
like the case of TR, the photon emission produced by 
a superluminal optical perturbation is not limited by 
the maximum variation of refractive index produced by 
optical perturbation. Instead, when the phase velocity 
of the waves <f>i and <j) r along u are identical to it, 
corresponding respectively to (3rii cos Q\ es — * 1 and/or 
(in r cos#£ es — * 1, the optical perturbation and the waves 
4>i and (f> r move together and can interact for longer times 
producing an arbitrarily large number of photons. This 
process can be described as a form of superluminal res- 
onance. We identify three regimes: i) for f3rii < 1 and 
(3n r < 1 there are no resonances; ii) for either (3rii > 1 
or j3n r > 1 there is only one pair of resonant emission 
angles; and iii) for both f3rii > 1 or (3n r > 1 two pairs 
of resonant light are emitted. These resonances can be 
achieved for a wide range of experimental parameters and 
configurations. The angular distribution corresponding 
to equation (|lip is represented in Figure 2 where we can 



Finally, we briefly discuss entanglement degradation 
due to thermal photons by considering that the system 
is initially in a equilibrium state with an average number 
of photons higher than the vacuum. In this case, the 
initial state is no longer the vacuum but rather a Gibbs 
thermal state of both initial wave modes, <pi and <f> a , i.e. 
Pia(n) = p(n) p(n), with 



P(R) = JL W_!L 
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where n is the thermal mean occupancy. After the optical 
perturbation has passed, the state ptriz) describing the 
4>t and <j) r modes reads 



p tr {z) = S(z)p la (n)S(z) 



(15) 



which is a squeezed thermal state [26J and for which Ey n 
is not an adequate entanglement measure [2?| . How- 
ever p tr is a Gaussian state, and its entanglement can 
be completely characterized using continuous variable 
methods (see 28] for a review), namely via the loga- 
rithmic negativity, Ejy(p) — max [0, — In p] , where p is 
the smallest sympletic eigenvalue of the Gaussian state 
p tr . The expression for p (see 29] for a derivation) is 
p = (2n + 1) exp (— 2arctanhz). The latter defines a ther- 
mal occupancy n c , above which all entanglement van- 
ishes, yielding 



2n c + 1 = exp (2arctanhz) . 



(16) 



Close to resonance (z — ► 1) the maximum allowable ther- 
mal occupancy diverges n c — > oo; entailing that entangle- 
ment extraction from optical boundaries is very robust 
regarding temperature by choosing a sufficiently high 
squeezing parameter. 
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In conclusion, we have presented a new emission mech- 
anism of entangled radiation using a sharp optical pertur- 
bation with an apparent superluminal velocity. The emis- 
sion spectrum and the emissivity depend on the apparent 
velocity and the change of refractive index of the opti- 
cal perturbation. These results extend those of reference 
[l3| from a one dimensional configuration to include all 
complex plane wave modes in a three dimensional space 
and valid for an arbitrary dispersive medium. For our 
particular configuration, the optimum direction of emis- 
sion is defined by the resonances /3rij cos Q\ es — > 1 and 
(3n r cos Q r r es — ► 1. The resonance angles Q r t es and d r r es cor- 
respond to both the best radiance and to the optimally 
entangled photons. From a purely theoretical point of 
view, this process has considerable advantages over PDC 
as a source of entangled light, namely since it is capa- 
ble of delivering two well separable and highly entangled 
beams with large intensities. In our case the photons are 
entangled in momentum whereas in PDC the photons 
are entangled in polarization, however these two types of 
entanglement can be interconverted [30l |. From a more 



experimental point of view, it is not easy to produce a 
sharp and sudden optical perturbation at scales inferior 
to the optical wavelengths to allow the large number of 
photon pairs necessary to make this process competitive 
with PDC. Present day technology allows to produce ion- 
ization fronts with 7 factors up to the order of 10 ± 10% 
that behave as a highly reflecting sharp optical boundary 
capable of exciting ZPF in the THz frequency domain and 
emit in the optical range [3l[. Moreover, a recent work 
has shown that this quantum mechanism of extracting 
photon pairs out of ZPF can be extended to optical per- 
turbations with arbitrary shape as long as they have an 
apparent superluminal velocity [Hj]. These results sug- 
gest the possibility of generating entangled photons via 
dynamical Casimir effect and can be generalized to arbi- 
trary time-dependent optical perturbations. 
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